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Abstract. We prove that the evolution of marginals associated to the empirical measure of 
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1 Introduction 

The hard-sphere dynamics plays a central role in the theory of dynamical systems and 
as underlying microscopic model in kinetic theory. The aim of this paper is to study 
some features of the time evolution associated to the finite hard-sphere model, namely we 
consider N identical hard spheres of diameter e moving in the whole space H-S system 
in the sequel. Let 

z = , 2 :* = , i = l,2, •■■A, 

be the initial configuration of the system. Here we denote by {xi,Vi) the position and 
velocity of the i-th particle. The evolution is deterministic and given by the usual laws of 
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elastic reflection. Let 


z(t) = , Zi{t) = {Xi{t),Vi{t)) e 




be the configuration at time t. On a full-measure set of initial configurations, the flow 
z —)■ z(f) exists for all times pQ. 

Choosing an initial conhguration z for which the dynamics is well posed, one dehnes 
the empirical measure as 


1 


N 


/iAr(t, dz) = — '^6{z - Zi{t)) dz , 


( 1 . 1 ) 


i=l 


where ^(■) is the Dirac measure at the origin (we drop the dependence of /iat on z). The 
sum of Dirac masses is transported through the H-S flow. Observe that the evolution of 


(1.1) is simpler than the evolution of, say, measures, since it amounts to follow one 


single trajectory. 

Is there a partial differential equation describing the evolution of /itv? Of course the 
empirical measure will not satisfy any closed equation, but rather a hierarchy of equations 


involving the higher order “empirical marginals” associated to (1.1), exactly as in the 


case of absolutely continuous distributions which evolve through the hard-sphere BBGKY 
hierarchy. The latter hierarchy has been studied in detail by several authors [SlIHKITlEoKISl 
[7j. In contrast with the well-known BBGKY holding for smoothly interacting systems, its 
derivation is delicate because of the singular character of the H-S flow. Indeed the collision 
operators appearing in the hierarchy are dehned by integrals on manifolds of codimension 
one, while the H-S flow is dehned only away from a null-measure set. This problem is even 
more delicate when one tries to give a meaning to the hierarchy for singular measures of 


type (1.1), since the integration of delta functions is performed on the boundary of the 


space where the H-S how takes place. 

In this paper, we show how to establish and rigorously justify the H-S hierarchy for em¬ 
pirical marginals. To do this, we make use of a series solution representation as introduced 
hrst in ra. where any integration over boundaries is carefully avoided. In particular, we 
show that among all the terms appearing in such representation, the “real” H-S trajectory 
is singled out by a mechanism of cancellations. 

As in the case of non-singular measures, the result can be obtained with several diherent 
methods. We point out that the method of reference (TB] does not require any regularity 
property and therefore can be easily extended to a general measure including singular 
parts. Approximation with smooth measures can be also used, together with some detailed 
information on the H-S flow. In the present paper we shall use a different approach. 
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which is simpler and more natural in the case of empirical marginals. The ingredients of 
the proof are the semigroup property of the flow and the boundedness of the number of 
collisions. Moreover, once established the result for empirical marginals, the validity of the 
H-S BBGKY hierarchy for absolutely continuous measures of class can be recovered by 
computing expectation values. Hence, the method of this paper can be also seen as an 
alternative approach to the validation issue discussed in the quoted references. 

Connected with the above analysis is the problem of describing the so called microscopic 
solutions to the Enskog equation, on which we comment in the last part of the paper. In the 
literature, several versions of the Enskog equation may be found. Here we will consider the 
form sometimes called Boltzmann-Enskog equation. This is a kinetic equation in which, 
in contrast with the Boltzmann equation, the diameter £ of the particles enters in the 
expression of the collision operator. More precisely, it reads as 


{dt + V ■'Vx)f{x,v,t) = X ^ dviduj {v — vi) ■ 00 (1.2) 

- Uje,v[,t)f{x,v',t) - f{x + , 

where the unknown / denotes the probability distribution of a given particle. As usual x, v 
and t denote position, velocity and time respectively. Moreover 5^ = {ca G S^\ {v — vi)-u > 
0}, S'^ is the unit sphere in (with surface measure du), {v,vi) is a pair of velocities in 
incoming collision configuration and {v', n() is the corresponding pair of outgoing velocities 
defined by the elastic reflection rules 


v' = V — oj[oj ■ {v — Ui)] 
v'l = Ui + oj[uj ■ {v — Ui)] 


(1.3) 


Finally A is the mean free path. Note that the Boltzmann equation is recovered, at least 
formally, for e = 0. In this paper, we will keep fixed £ > 0. 


In 1975 N.N. Bogolyubov |1] observed that there exist solutions to eq. (1.2) of the 
form (1.1), where {ziit)}^^^ denotes the evolution of an Y—particle H-S system. These are 


called microscopic solutions. 

Let US add here a comment concerning this terminology. We recall that, for N large 


and corresponding small e, the (1.1) should provide the normalized density of particles in 
the microscopic description of the gas evolving along (1.2). Considering z as a random 


variable and under a suitable “chaotic” assumption, (1.1) can be actually shown to be very 

However the convergence 


close to the solution of the Boltzmann-Enskog equation 
takes place only in the continuum limit where Y —)■ cx) as 
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Bogolyubov’s statement concerns instead (1.1) and (1.2) for finite N and £. The re¬ 
sult may look surprising, the Enskog equation being a genuine irreversible kinetic equa¬ 
tion. To see this, one may introduce the free-energy functional H{f) = f flogfdxdv + 
/ir 3 Ib(x e) d(^) Phj) ) where p is the spatial density and i?(x, e) is the ball around 
X and of radius £. Then it turns out that H decreases if / is a solution to (1.2), see e.g. 
[3]. But since the functional H does not make sense on solutions of the type (1.1), there is 
no a priori contradiction. 

As in the case of the H-S hierarchy, the Boltzmann-Enskog collision operator appearing 
in the right hand side of (1.2) is not well dehned when evaluated in / = pM, so that 
a discussion on the precise mathematical meaning of the microscopic solutions of [1] is 
required. In [THl [IS] a suitable regularization of has been used to give a sense to (1.2) 


in terms of a limiting procedure. In this paper, we will approach the problem in a different 
way. 

Motivated by the fact that the concept of series solution appears to be convenient to 
justify the microscopic solutions to the H-S hierarchy, we shall focus on the same notion of 
solution for the Boltzmann-Enskog equation and on the comparison between them. We hnd 
it interesting to observe that such a comparison is rather non-trivial. The series solution of 
the H-S system provides a unique well dehned result, when applied to microscopic initial 
data. In contrast, the corresponding expansion for the Boltzmann-Enskog equation does 
not allow to avoid the integration of singular measures over boundaries. As a consequence, 
not only the single term of the expansion depends on the regularization, but the series 
is not even absolutely convergent for short times. The origin of the ambiguity is that 
the space of singular solutions of the kinetic hierarchy allows contractions, i.e. different 
particles having the same conhguration. Unfortunately, a natural prescription preventing 
this phenomenon seems to be missing. 


The plan of the paper is the following. In the next section we introduce the H-S and 
Enskog hierarchies for data, together with the tree expansion, which is a basic tool 
for our analysis. In Section we introduce the microscopic states, namely states which 
are concentrated on a single conhguration, and extend the notion of hierarchy to this 
context. After that, we analyze the series solution to the Boltzmann-Enskog equation for 
microscopic initial data. 


4 









2 Hierarchies 


In Section 2.1 below we introduce the H-S system, recall the preliminary results on the 
H-S dynamics and explain how to describe the evolution of a class of absolutely continuous 
measures, by means of a hierarchy of equations similar to the usual BBKGY hierarchy 
for smooth potentials. An analogous description is also given for the Boltzmann-Enskog 


evolution. In both cases, we provide in Section 2.2 the explicit representation of the series 
solution in terms of the tree expansion and of a class of special flows of particles evolving 


backwards in time. We mainly follow Sec. 2 of ref. na in this part. Finally in Section [273 


we show how to properly formulate the series solution in order to extend it to the case of 
singular measures. This requires a discussion on the invertibility properties of the flows. 


2.1 Preliminary results 

We consider a system of N hard spheres of unit mass and of diameter £ > 0 moving in 
With Zi = {xi,Vi) G we indicate the state of the iWh particle, i = 1, 2, ■ ■ ■, while for 
groups of particles we use the notation Zj = (^i, • • • , Zj) . “Particle i” is a particle whose 
conhguration is labelled by the index i. 

The particle conhguration lives in the phase space of the system, dehned as the subset 
of in which any pair of particles cannot overlap, namely 


M-n — ^ G 


d6A^ 


\Xi 


Xfcl > 0, = k jk 


( 2 . 1 . 1 ) 


Given a time-zero conhguration zjv G Ad at, we introduce the how of the Y-particle 
dynamics (equations of motion) 

1 1—)■ T7v(t) zn , (2-1.2) 

by means of the following prescription. Between collisions each particle moves on a straight 
line with constant velocity. When two hard spheres collide with positions Xj, Xj at distance 
£, normalized relative distance 


uj = {xj — Xi)/\xi — Xj\ = {xi — Xj)/e G S‘^ 

and incoming velocities Vi, Vj (i.e. {vi — vj) ■ a; < 0), these are instantaneously transformed 
to outgoing velocities n',n' (i.e. {vi — Vj) ■ a; > 0) through the relations 

v[ = Vi - u[oj ■ {Vi - Vj)] , 

n' = Vj + u[u ■ {vi - Vj)] . (2.1.3) 
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Such a collision transformation is invertible and preserves the Lebesgue measure on M®. 
Notice also that the flow T 7 v(t) is piecewise continuous in t 

The above prescription does not cover all possible situations, e.g. triple collisions (three 
or more particles simultaneously at distance e) and grazing collisions ((u* — Vj) - a; = 0) are 
excluded. However what is important is that the flow is globally dehned, almost everywhere 
in In fact we have from [1] (see also [II1|6]) that there exists in Aljv a subset, whose 
complement is a Lebesgue null set, such that, for any z^v in the subset, the mapping (2.1.2) 
is a solution of the equations of motion having TAr(O) z^v = Moreover, the shifts along 
trajectories t ^ T^it) dehne a one-parameter group of Borel maps on AItv which leave 
the Lebesgue measure invariant. 

In particular, it follows that the number of collisions in any hnite interval is hnite for 
almost all initial conhgurations. Actually, in our situation one has a stronger result mm-- 


Proposition 1 There exists N E N such that the total number of collisions in the H-S 
flow is at most N. 

We shall make use of this property in the present paper. It is worth stressing that everything 
that follows would automatically apply also to N spheres enclosed in a region A C with 
elastically reflecting boundaries, as in [Tj, in which case the above proposition is generally 
true in any hnite time interval only after removing a suitable Lebesgue null set (one single 
particle in a box may undergo inhnitely many collisions in a hnite time). 

A statistical description of the system is provided by assigning on Ad at an absolutely 
continuous probability measure with density W, initially - and hence at any positive time - 
symmetric in the exchange of the particles. Its time evolution is described by the Liouville 
equation, which in integral form reads 


lV(zjY,t) — lVo(TjY(—t) Zjy) , 


(2.1.4) 


a.e. in A4jy, where Wq is the assigned initial datum. The j-particle marginals for j = 
1, ■ ■ ■ ,N are given by 


fj{zj,t)= dzj+i-■ ■dzNW{zj,Zj+i,-■ ■ ,ZN,t) , (2.1.5) 




where Zj E M. j and 


S{'Kj) = < z = {x,v) E \x — Xk\ > e for all/c = 1, • • • , jA . (2.1.6) 


6 




Moreover, to simplify the following notations we shall extend the dehnition over by 


= 0 , 


Zj e \ Mj. 


(2.1.7) 


The evolution equations for the considered quantities were hrst derived formally by 
Cercignani in and take the form (II-S BBGKY hierarchy) 


{dt + Cj) fj{zj,t) = / duj dvj+i B{uj-,Vj+i-Vi) fj+i{zj,Xi + euj,Vj+i,t) , 

,_i 2s2xr3 

( 2 . 1 . 8 ) 


where £,■ denotes the generator of the j—particle dynamics dehned by 


snt)fAz„-) = 




Zj G Mj 


G \ Mj 


(2.1.9) 


and the collision kernel is 


B(u; Vj+i - Vi) =u- (vj+i - Vi) . 


( 2 . 1 . 10 ) 


Notice that, by virtue of (2.1.7), the integration du is restricted over the subset of the 
sphere 

{a; G I min \xi + cue — xi] > e} . (2.1.11) 

£=!,■■■ ,j;£^i 

The series solution of the hierarchy, obtained perturbing the j—particle evolution, is 

rt rh rtn-i 

IQ JO 


- j,n) dti dt2-- - I dtn 
n=0 Jo Jo Jo 

Sj{t — tl)Cj+lSj_^_-^^{tl — ^ 2 ) ■ ■ ■ Cj+n<Sj_^_nitn)fo,j+n , ( 2 . 1 . 12 ) 

where /j(-,0) = /oy is the initial datum (marginals of Wq), 

a(r, n) = r(r — 1)... (r — n + 1) (2.1.13) 

{a{r, 0) = 1), and Cj+i is the collision operator, given by the sum in the right hand side of 


Eq. (2.1.8^, i.e. 




*J+1 


i=l 


CiJ+lfj+l {Zj , t) — 


'52xR3 


(2.1.14) 

du dvj+i B{u-, Vj+i - Vi) fj+i{zj, Xi + eu, vj+i, t) . 
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The rigorous validation of formula (2.1.12) has been discussed in [IZIEIESIIZI and can be 
proved under rather weak assumptions on the absolutely continuous initial measure. An 
alternative proof will be presented in Section]^ (Corollary [^. 

Finally, we introduce the so called Enskog hierarchy. 


Let 5 ^ be a solution to the Boltzmann-Enskog Equation (1.2). Then the products 


satisfy 


:= g{t)®^{7.j) 


dt T ^ ^ 'Vi ■ ^Xi I 9j ^ ^j+^9j+i • 

i=l 


(2.1.15) 


(2.1.16) 


This can be easily obtained performing a change of variables u —>■ —u inside the positive 


part of the collision operator in (1.2). The corresponding series solution is: 

fii rtn-i 

/o 


djit) = X] / dt2-- - I dtn 

n>0 do Jo 

•Sj[t ti)Cj-|-nSjf-|_i(ti ^ 2 ) ■ ■ ■ ^j+n^j+n(tn')90,j+n ; 
where now Sj{t) is the free flow operator, dehned as 

‘5i(^)/i(Zj, ■) = fj{Xi - UiLUi, • • • ,Xj- Vjt.Vj, •) , 

and 


(2.1.17) 


(2.1.18) 


(2.1.19) 


9oj = 9o 

is the family of initial data. Existence and uniqueness for the solutions to the Enskog 
hierarchy have been discussed in |2l |3l |13] . 


In spite of their formal similarity, the two expansions (2.1.12) and (2.1.17) are deeply 
different as explained in [T3] and as conhrmed also by the discussion in Section 3.1 below. 


Let US conclude here by establishing a fundamental property of both the series expan¬ 
sions introduced above, namely the semigroup property which is a consequence of the same 


property holding for the operator (2.1.9) and can be formulated in the following way. 


First we introduce the operator 7^(zj,t) acting on the marginal of order j + n and 


describing the n—th term of the expansion (2.1.12): 


rt pti rt„-i 

Tnizj, t)fj+n = aiN - i, n) dti dt 2 --- dtnS^{t - ti)Cj+i ■ ■ ■ . 

Jo Jo Jo 

( 2 . 1 . 20 ) 










For reasons that will be clear in the next section, we may call Tn{zj,t) the n—particle tree 
operator for the interacting flow. Eq. (2.1.12) can be written 


N-j 

= '^Tn{Zj,t)fj+n{0) . 
n=0 

If t + t 2 , by simple algebraic manipulations it follows that 

N-j 

= '^Tnizpti)fj+nit2)- 

n=0 

In a similar way for the Enskog flow we define the n—particle tree operator 

rt rti rt„-i 

7^ (Zj, t')Pj^n / dti I dt 2 ■ ■ ■ / dt^Sj(t ti')Cjj^i • • • 9j+n 

Jo Jo Jo 

so that 

oo 

9j{2j,t) = '^rj^{zj,t)gj+n{0) 

n=0 

and there holds 

OO 

n=0 


2.2 Tree expansion 


In this section we shall write formulas (2.1.12) and (2.1.17) in a convenient 
explicit way. We follow [151 E]. 


Extracting the sums from (2.1.14) in formula (2.1.12), one has 

f*^n —1 


V—^ ^ I /■* /■*! 

/i W = / ^^2 ■ ■ ■ 


dtr. 


n=0 k„ 


where 


■<Sj{t t2)---Cknj+n<Sjj_n{tn)fo,j+n, 


j i+1 j+n-1 

k„ ki=l k2=l kn = l 


To describe the summation rule we introduce a useful notation. The n- 
j—particle tree, denoted by r(j, n), is defined as the collection of integers ki, ■ ■ 
are present in the sum (2.2.2[), i.e. 


ki G Ij, k2 G dj+ii • • • ,kn^ Ij+n-l , 


with 


4 = ,4- 


( 2 . 1 . 21 ) 

( 2 . 1 . 22 ) 

, (2.1.23) 

(2.1.24) 

(2.1.25) 

and more 

( 2 . 2 . 1 ) 

( 2 . 2 . 2 ) 

-collision, 
■ , kn that 

(2.2.3) 
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In this way 


(2.2.4) 


k„ r(yn) 

The name “tree” is justified by its natural graphical representation, which we explain 
by means of an example: see Figure corresponding to r(2, 5) given by 1,2,1, 3,2. In the 



Figure 1: The two-particle tree F(2,5) = 1, 2,1, 3,2. 

figure, we have also drawn a time arrow in order to associate times to the nodes of the 
trees: at time the line j + i is “created”. Lines 1 and 2 of the example, existing for all 
times, are called “root lines”. 

Note that a j—particle tree is not a collection of j one-particle trees because, in the 
latter case, the ordering of particles belonging to different one-particle trees is not specified. 
Given a j-particle tree F(j, n) and fixed a value of all the integration variables in the 


expansion (2.2.1) (times, unit vectors, velocities), we associate to it a special trajectory 
of particles, which we call interacting backward flow (IBF in the following), since it will 
be naturally defined by going backwards in time. The rules for the construction of this 
evolution are explained as follows. 

First, we introduce a notation for the configuration of particles in the IBF, by making 
use of Greek alphabet i.e. where s G [0,f] is the tim^ Note that there is no 

label specifying the number of particles. This number depends indeed on the time. If 
s G {tr+i,tr) (with the convention = t,tn+i = 0), there are exactly j + r particles: 


C('5) (Cl(®)) ' ' ' )C;+r('S)) G .Arlj^r foi’ S ^ (G+IjG) ) 


(2.2.5) 


Bn previous work mun] the notation C^{s) has been used for the IBF (and C(s) for the corresponding 
“Boltzmann flow”). Here we drop the superscript, since we will keep e fixed throughout the paper. 
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with 

Ci(s) = te(s),^i(s)) , (2.2.6) 

the positions and velocities of the particles being respectively 

^(s) = (6(s),--- , W(s)) , 

r]{s) = (r/i(s),..., Vj+r{s)) • (2.2.7) 


Our hnal goal is to write Eq. (2.2.1) in terms of the IBF (to be dehned below). More 


precisely, (2.1.20) shall be rewritten as 


'T~niCji^)fj+n <y{N j,Tl) ^ ^ / dA(tn, U^ni ^j,n} J_ J_ (tj))/j+n (C (0)) ) 

rij,n) i=l 

( 2 . 2 . 8 ) 

where = C,j{t) := , (C(s))sg[op is dehned below, dL is the measure on M"' x x 


dA(tjj, ttJnj j,n) • • • dtjiduJi . . . duJ■)ldVj-^-l . . . dVjjj-n , 

with the abbreviation 

(fnj ^71) ^j,n) (^1) ^25 ■ ■ ■ 5 ^1; ' ' ' 5 0+1’ ' ' ’ ’ 0+^) ’ 


(2.2.9) 


while 


B [uJi^Vjj^i Vkiiti)) B[uJi,Vj^i hfci(^*))ii{|^j+i(ti)-?fe(ti)|>£ (2.2.10) 


with B dehned by (2.1.10). Summarizing, the term Tn{C,ji ^)/j+n(0) is the sum over all trees 
r(j, n), of terms where the initial datum /o,j+n is integrated, with the suitable weight, over 
all the possible time-zero states of the IBF associated to r(j, n). 

In formula (2.2.8), the triple (f,, cuj, Uj+j) is thought as associated to the node of r(j, n) 
where line j+i is created (see Figure[^. In the rest of the paper, we shall further abbreviate 


f dA{tn,U}n,^j,n)YlB^ = / dA(t„, V^'„) JJ 5' 

'' i=l 


{uJi,Vj^i 'Qkiid'i)') ) (2.2.11) 


where the rjkX'ti) in factors B have to be computed through the rules specihed below, 
starting from the set of variables (t„, Vj_„), the corresponding j-particle tree (whose 

nodes are labeled by (t„,Vj^„)), together with the associated value of Cj,t. 

Let US hnally construct <^(s) for a hxed collection of variables T{j,n), Q,tn,uJn,^j,n, 
with 

t = to > ti > t2 > ■ ■ ■ > tn > tn+1 = 0 , (2.2.12) 
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and ujn satisfying a farther constraint that will be specihed soon. The root lines of the 
j-particle tree are associated to the hrst j particles, with conhgnration Cir'' yCj- Each 
branch j + i {i = 1, ■ ■ ■ ,n) represents a new particle with the same label, and state (j+e- 
This new particle appears, going backwards in time, at time ti in a collision state with a 
previons particle (branch) G {1, • • ■ j + £ — 1}, with either incoming or ontgoing velocity. 

More precisely, in the time interval {tr,tr-i) particles 1, • ■ ■ , j + r — 1 flow according 
to the nsnal dynamics This dehnes starting from Cj+r-iitr-i)- At time tr 

the particle j' + r is “created” by particle kr in the position 




(2.2.13) 


and with velocity Vj+r- This dehnes C(A) = (Ci(A), • • • , Cj+r(tr))- After that, the evolntion 
in {tr+i,tr) is contrncted applying to this conhgnration the dynamics Tj+r (with negative 
times). 


The characteristic fnnction in (2.2.10), is a constraint on ojr ensnring that two hard 


spheres cannot be at distance smaller than e. 

We have two cases. If ujr ■ {vjj^r — Vkritr)) < 0, then the velocities are incoming and no 
scattering occnrs, namely after tr the pair of particles moves backwards freely with veloc¬ 
ities rjk^itr) and fj+r- If instead Ur ■ {vj+r — Vkritr)) > 0, the pair is post-collisional. Then 
the presence of the interaction in the how Tj+r forces the pair to perform a (backwards) 
instantaneous collision. The two situations are depicted in Fig. 



Figure 2: At time U, particle j -|-r is created by particle kr, either in incoming {oor ■ {vj+r — 
Vkritr)) < 0) or in outgoing (cn^ • ivj+r — Vkritr)) > 0) collision conhgnration. Particle kr is 
called “progenitor” of particle j + r. 

Proceeding inductively, the IBF is thus constructed for all times s G [0,f]. 

Remark 1 Between two creation times tr,tr+i any pair of particles among the existing 
j + r, diherent from the couple (fcr, J + ?’), can possibly interact. These interactions are 
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called recollisions, because they may involve particles that have already interacted at some 
creation time (in the future) with another particle of the IBF. In our language, recollisions 
are the “interactions different from creations”. 

To conclude, we have obtained the following representation for the series solution to 
the HS-BBGKY hierarchy: 

N-j 

/i(ziG) = a{N -j,n) 

n=0 r(j» 

For future purpose, we write also the version of the expansion where we split positive 
and negative part of B by setting 


(iA(t 


ni ^ri’) ^ 




(2.2.14) 


= B% + Bi_ 


with 

Bl{uj,V) = B^{uj,V)U^.v>o . 


If (Tn = (cTi, ■ ■ ■ , an) where cxj = ±, one has 


N-j „ 

= X] X] - j,n) j dA(t„,n;„,Vj-„) Jj5^^/oj+„(C(0)) , (2.2.15) 


n=0 r(i,n) 


where = Y\a=i ~ Vkiiti)). We underline that, according to our previous 

definitions, Bt gives a negative contribution to the series expansion (2.2.15). 

We also notice that the trajectories entering into expansion (2.2.15) have nothing to 
do, in principle, with the real trajectories performed by the particle system. However the 
latter can be recovered by the expansion, by means of a complex system of cancellations. 
The forthcoming discussion on the hierarchy for empirical measures will clarify this point. 


Finally, a similar analysis can be done for the Enskog hierarchy. We introduce the 
backwards flow C^{s), called Enskog backwards flow (EBF), which is constructed as <^(s) 
with the additional prescription that all the recollisions are ignored. In this flow, particles 
are still created at distance e from their progenitor, but they may overlap, i.e. they may 
reach a distance smaller than £ during the evolution. In particular, the time-zero state 
<^®(0) lies in 

Alternatively, we may say that the EBF is constructed exactly as the IBF, except for 
the following differences: 

- the interacting dynamics T is replaced by the simple free dynamics; 

- there is no constraint on ujr. 
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In terms of this flow, Eq. (2.1.17) can be written explicitly 


OO « 

J^i?o-„5'0J+n(C (0)) • (2.2.16) 

n=0 r(j,n) <Tn 

where H = nr=i ^j+i “ VkS'^i)) and V) = B{u;, V)l±^.v>o ■ 

It is not difficult to check that the EBF allows a complete factorization, whenever 
the initial datum does na. Namely if is taken as in ( |2.1.19[ ), then ( |2.2.16[ ) gives 

gj{zj,t) = gi{t)^^{zj). 


We end this section by summarizing the differences between (2.2.15) and (2.2.16): 

• vs. E 


^n>0 

a{N — j, n) vs. A” 
Bt vs. B„ 

(Tn 

fo,j+n VS. 5^0,j+n 

IBF vs. EBF 


(total number of particles of the system); 
(multiplicative coefficient in the expansion); 
(overlap of created spheres prevented / allowed); 
(initial data, respectively non-tensorized / tensorized); 

(backwards flow with recollisions / overlaps). 


In particular, we observe once more that, in contrast with g^ j which can be taken as a 
product state, foj cannot factorize because it must prevent the overlap between different 
spheres (condition (2.1.7)). 


2.3 Weak formulation of hierarchical expansions 


In this section we discuss a weak formulation of the above introduced series solution to the 
HS hierarchy, suitable to deal with singular measures. 

Let us consider the right hand side of (2.2.15) in the case that the initial data fo,j+n 
are replaced by suitable measures, not absolutely continuous. To give a precise meaning to 
the formula, we need to integrate with respect to the initial conhguration variables C(0). 
This amounts to study in some more detail the IBF-map dehned, for hxed f > 0, r{j,n) 
and (Tn, by 

Cji^n, t ^j-|-n,(0) , (2.3.1) 


where Cj = zj. It follows from the properties of the H-S flow that the above transformation 
is a Borel map almost everywhere dehned over the domain specihed in the previous section, 
with image A{t, T, cr„) C see e.g. Lemma 1 in ng , or 1201 for a different method 

^I.e., dCj+n—essentially, 

A(t,r(j,o),0) = M,+„, 
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of proof. Moreover, a simple computation shows that the Jacobian determinant of the 
transformation is in modulus nr.i Vj^i — Tjki{ti))\, SO that the map induces the 

equivalence of measnres 




(2.3.2) 


We constrnct now the inverse of the map. To do so, we introdnce the interacting forward 
flow (IFF) 

Ci+n(0) )■ (C ('S))^gjQj] 

dehned in the following way. As for the IBF, the IFF has a nnmber of particles which 
depends on time. At s = 0, one has j + n particles. Let us take for such particles a 
conhguration <^j+„(0) in the image A of (2.3.1) (which depends on t, T{j,n) and cr„). We 
let the conhguration evolve forward in time via the H-S dynamics Tj+n up to the hrst 
collision between particle j + n and particle {j + n is the last one created in the IBF 
and generated by particle kn according to the tree T{j,n)). Such an instant of time exists 
because <^j+n(0) belongs to A. There are two possibilities. 


• This interaction is a creation. Then particle j + n disappears, while particle kn 
interacts or not according to = + or cr„ = — respectively. Next, particle 
evolves together with the other j + n — 2 particles in the H-S dynamics Tj+„_i, np to 
the next collision dictated by F(j, n), namely the collision between particles j + n — l 
and kn-i- 


• This interaction is a recollision. Then we let evolve the system fnrther with the 
dynamics np to the next contact between j + n and kn- Clearly, meanwhile 
other recollisions occnr. If there is no next contact in [0, t] between j + n and k^ 
we simply eliminate this second option. 


As we shall see in Section 2.3.1 there are cases in which both options are possible. There¬ 
fore, the iteration of the above procedure generates M different flows i = 1, - , M, 

with M depending on <^j+„(0). Of course by Proposition]^ there exists W G N such that 
1 < M < 2^. 


We conclude that (even though locally invertible) the map (2.3.1) is not globally in¬ 
vertible and integrating Eq. ( |2.2.15 ) against a bonnded continuous function (f = ■ 

A(t,r(j, 1) = fci,cri) = |o+n+i(0) G Mj+n+i “j + n + 1 and ki collide in (0,t)”| 

and so on. 

^More generally, if iteration of the procedure becomes impossible. 
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—>■ M, the following weak formulation of the time-integrated H-S BBGKY hierarchy is 
obtained: 


„ N-j n 

j= E Eli 

n=0 T{j,n) cr„ r=l 


M 


ar a{N — j, n) e 


—2n 


' A{t,r,cr„) 


dCj+n fo,j+n{Cj+n) E (*)). 


i=\ 


(2.3.3) 

where M = : A —>■ N. Note that the latter function is an a.e.-dehned step function 

over Aij+n, strongly dependent on the details of the H-S flow. Observe also that the right 
hand side of (2.3.3) makes sense even for initial marginals which are measures d/oj(Ci)- 
Though the above formula is not very handable for practical purposes, we shall use it 
only in the case of discrete measures. Remind that by (2.3.2), the transformation (2.3.1) 
is nonsingular out of the grazing collisions = 0. Therefore for any Dirac measure 
supported in points of H(f, T, cr„) “avoiding” grazing collisions, (2.3.3) can be rewritten 
without ambiguity as 


N—j 

J /.,(*¥ = E E E“(^-J'") / n 


Wj-Ate). (2.3.4) 


n=0 T{j,n) cr„ 


In other words, we describe the discrete measure as a pushforward measure along the 
mapping ( |2.3.1[ ). 

In Section]^ we will show that the empirical marginals associated to (1.1) and supported 
on “good” conhgurations of the H-S flow satisfy indeed (2.3.3). We shall call these solutions 
the microscopic solutions of the H-S hierarchy. An obvious way to proceed would be by 
taking a sequence of data (/oj)fc approximating the singular measure as fc —)■ cx) and using 
some topological information on the function M = M(<^j+„). Instead, we shall present in 
Section a more constructive approach. 

We conclude by observing that the EBF-map dehned, for hxed f > 0, r(j, n) and cr„, 
by 






Cf+n(0) 


(2.3.5) 


with = Zj, is globally invertible over its domain of dehnition. This follows immediately 
from the fact that each pair of particles in the EBF may interact at most once. As a 
consequence, from (2.2.16) we easily deduce the following weak formulation of the time- 
integrated Enskog hierarchy. 


P oo n 

1 9j(i¥ = E E Eli 

n=0 r{j,n) o-n r=l 




lAE(t,r,cr„) 


dCf+n 90,j+n{Cf+n) 0(Cf (^)) , (2-3.6) 
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where A^{t,r,crn) C is the image of (2.3.5). In Section 3.1 we show that this 


formula is not well-dehned for the empirical measure, since ambiguities arise from the 
contractions in goj+n and the corresponding integral of the singular measure over the 
boundary of A^. 


2.3.1 Non-invertibility of the IBF 


We show the non-invertibility of the map (2.3.1) by means of an example 


Let us consider the case of = 3 hard spheres, with initial conhguration C 3 ( 0 ) such 
that the dynamics admits the sequence of collisions as in the hgure that follows, i.e. Ci = 
(2,3),C2 = (1,2),C3 = (2,3),C4 = (1,2). 



(2.3.11) 


Many conhgurations of this type can be constructed in two or three dimensions [12]. The 
above hgure is an IBF for the following tree r(l, 2) with positive nodes <72 = (-I-, +): 



(2.3.12) 


The collisions ci and C 2 correspond to the nodes of the tree, while C 3 and C 4 are recollisions. 

A different IBF for the same F(l, 2) and <T 2 is obtained by taking C 3 and C 4 as creation 
times, as in the hgure below which, for the same initial conhguration and tree, yields a 
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(2.3.13) 


3 Microscopic solutions to the H-S hierarchy 

Let us start by introducing the probability measure over Mn having as 1—particle marginal 


the empirical measure. The density W and marginals fj introduced in Section 2.1 for the 
absolutely continuous case are here replaced by the empirical marginals, which are higher 
order probability measures concentrated on subsets of the H-S configuration. 

Let zjy = {zi,--- ,Z]si) G be a given configuration. A symmetric probability 

measure over Ad at concentrated on z^r is given by the empirical density 


N 




TT 2 = 1 


where vr is the generic permutation. 

The empirical marginals are dehned over Mj, as usual, by 


It follows that 


— J d(j+i ■ ■ ■ d(N ^(Cj5 Ci+i) • • ■ ) Cn) ■ 
Ai(Ci) = J d(2 - ■ ■ d(N A{(i, C 2 , • • • , Cn) 

N 

--J2 / dC2---dC^l[6{Q-z^(^) 

' IT d j=l 

^ ^ ^(Cl ~ ^7r(l)) 


1 

m 

1 

TV! 


N 


N 


Y.s{c^-z,), 


(3.1) 


(3.2) 


(3.3) 


2 = 1 
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namely 


^liOdC = HN{dC) 


is the empirical measure (1.1) at given time. By a similar computation: 


1 


^jiCj)— ^TTp)) 


m 

77 i=l 

{N-3)\ 


N\ 




j” ■ ^1 

1 


Af(Af - l)-- - (Af-j + 1) , 


E fl'SK.--*.) 


(3.4) 


^1 ?”■ ^1 


Observe that over Aij one has the identity: 


A.(C) = 




N{N-l)...{N-j + l) 




(3.5) 


Next we consider the time evolution of the marginals. 

First of all, we need to take care on how to fix the initial configuration zjy. By the 
results of roiniE] quoted in Section |2.1| one can identify the full Lebesgue-measure set of 
good H-S configurations 

M*^cMn (3.6) 

with the set for which triple (or multiple) collisions, grazing collisions and simultaneous 
collision^ are forbidden. The evolution of such configurations is defined for all times. The 
following set is of course still full: 


= {zw e Ml, zt e Ml Vz, c z„} . 


(3.7) 


From now on, we shall fix zjv € This is needed to ensure well dehned flows in the 

right hand side of (2.3.3)-(2.3.4). 

By dehnition. 


N 


A(Ck,«) = A(t„(-()Ck) = 


(3.8) 


77 2 = 1 


^Namely two pairs of particles colliding at the same time. This does not lead to an ill-defined dynamics, 
but is conveniently removed to simplify the following argument. 
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and (as in (3.4)) 




1 

N{N - 1)... {N - j + 1) 


2l,-” ,ij 5 = 1 


where here and below we indicate 


(3.9) 


Z7v(t) — TAr(t)zAr . 


We want to show that the seqnence formed by the empirical marginals of the time 
evolved measnre, satisfies the H-S hierarchy, Eq. (2.3.4). 


To do this, we nse Proposition and call S the nnmber of collisions delivered by ZAr(-) 
in the time interval [0,t). We partition [0,t) by a seqnence of S' + 1 intervals 

s 

[0,t) = lj[t„t,+i) (3.10) 

1=0 

where 0 = to < < ■ ■ ■ < ts+i = t and the following properties are satished: 

Property 1 In each time interval (tj, tj+i) with i = 0, - ■ ■ , S — 1, the H-S system delivers 
a single collision at time Tj+i. In the last interval [ts',t 5 +i) the motion is free. 

Property 2 If the collision at time Tj+i occnrs between particles a and /3, then tj+i — Tj+i 
is chosen so small that the following two different H-S flows are free: 

S -)■ TAr_i(s)(^^i(ti), • ■ • ■ ■ ■ , ZNiU )^ S G [0, ti+i - U ) 

and 

S -)■ TAr-l(s)^Zi(ti), • • • , Zfj - i { ti ), Zf ^+ i { ti ), ■ ■ ■ , ZN { ti)j s e [0,ti+i - U ) . 

Note that Property 2 means that if we remove particle a (or (]) jnst before the collision 
instant Tj+i, then we see a free flow np to tj+i. This will allow to restrict the compntation 


of (2.3.4) to n = 1,2. 


The existence of the above partition is gnaranteed by the assnmption zat G and 
by the continnity of the free flow. 


Now notice that, in view of the semigroup property Eq. (2.1.22), it is enongh to prove 


(2.3.4) for any time interval in the partition. Therefore in what follows we shall restrict 


to a time interval ( 0 ,f) in which at most one collision takes place and t is so small that 
Property 2 is fnlfilled. Withont loss of generality we set a = 1, /? = 2. 


20 








For simplicity of notation, we drop the integration f d(^j(j){Q) in many of the following 
formulas. Furthermore, we abbreviate Aj(0) = Aj. 

Let us consider hrst the case j = 1. 

The term n = 0 in the expansion (2.3.4) is 


ro(Ci,t)Ai= smMCi) 

1 ^ 

= {xi + Vit))6{r]i-Vi) . (3.11) 

i=l 

If S' = 0 (no collision is delivered by Z 7 v(-))) fhis reproduces trivially Ai(t). Otherwise 
if S' = 1 (one single collision occurs), all the terms in the above sum contribute to Ai(t) 
except those with i = 1, 2, which therefore must be cancelled by some other terms of the 
expansion with n > 0. 

If S' = 0, all such terms are easily shown to be zero. Let us assume S' = 1 and consider 
the term with n = 1 and ai = +, namely 


V(Ci,i)A2 

(iV-1) 


dA{ti,u, V 2 ) 5+ A2(C2(0)) 
[ dA{ti,u,V2) Bl 


(3.12) 


N{N-1) 

5(Ci(0) - Zi)5(C2(0) - ;^2) + ^(Ci(0) - ;^2)^(C2(0) 




Observe that only two of all the terms dehning A 2 appear in the right hand side of the 
equation, because all other pairs in zat do not interact by the above assumed properties, so 
that they fall away from the image of the IBF-map. In order to have other contributions it 
is then necessary that particle 1 (or particle 2), ignoring particle 2 (or particle 1), collides 
with some other particle. But such events are absent because of Property 2, invoked just 
to avoid these contributions. 


To compute the integral in (3.12) we proceed as in (2.3.3), namely we change variables 
according to Ci,ti, 0 Ji,V 2 —)■ Ci(0),C2(0) and use (2.3.2): 


<iCi^(Ci)7:+(Ci,«)A2 = ^ 


d(i d(2 ^(Cf(«)) i5(Ci - =i)i(C 2 - 22 ) + 5(0 - J2)5(C2 - 21 ) 


(3.13) 

where the integral on the right hand side is now extended over 24(t, F(l, 1),+), i.e. on 
the set of couples of particles leading to a collision. Since {zi,Z 2 ) lies inside this set, the 
integral is equal to 0 ( 2:1 (t)) + 0 ( 2 : 2 (t)). Therefore we obtain 


r+(Ci,t)A2 = - 


1 ~ ^i(^)) + <^(Ci ~ ^2 


(3.14) 
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With the same computation we get for the term with n = 1 and ai = —: 


Ti (Ci,^)A2 = - — <{5(6 - {xi + vit))6{r]i -vi) + 5(6 - {x2 + V2t))6{T]i - 


The term (3.15), which is negative, cancels the two terms in (3.11) with i 
we conclude that 


^ 2 )} . (3.15) 
= 1,2. Hence 


—To{Ci,t)Ai + T[^{(i,t)A2 + Ti {(i,t)A2 . (3.16) 


It remains to show that all the other terms in the expansion (2.3.3) with n > 1 vanish. 
This follows immediately from Property 1 and Property 2. In fact, in order to hnd a subset 
of 1 + n particles of zjv lying inside the image of the IBF-flow with n > 1, we would need 
at least one more collision with respect to those allowed by the properties. 

For j > 1, we proceed in the same way. The term To{Cj,t)^j is given by 


‘51(6A,(6) 


1 

N{N - 1)... {N - j + 1) 


E n'^(C-+.(^))- 


(3.17) 


As before if S' = 1 and particles 1 and 2 collide, then the terms which contribute to Aj{t) 
are those with either 


(1,2) c 61 ■■• 6 } 


or 

( i , 2 ) n 6 i ---*6 = 0 - 

The other terms, namely those for which 1 G {6 ■ ■ ■ +•} and 2 ^ {6 ■ ■ • or the reverse 
situation, are exactly compensated by 7^~(6', t)Aj+i, while the contributions missing in 
to reconstruct Aj{t) are produced by 7^^(6,t)Aj_|_i. Finally the same arguments 
used for j = 1 show that all the terms with n > 1 are vanishing. 

In conclusion, we have proved the following result. 



Theorem 1 Let zjv G The empirical marginals Aj{Q,t) defined in (3.4) and sup¬ 

ported in ztv verify, for any t > 0, the BBGKY expansion (2.3.3). Namely, for any bounded 
continuous (j) : —)• M, 
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By virtue of (3.5), this result can be formulated, in terms of the empirical measure /tat 


over A4j, as 


N-j 


(mwr = E E E 

n=0 r(j,n) (T„ 

where we set /iv( 0 ) = jujv, or, equivalently: 

„ N-j n 

I (M*)r> = EEEn 


2\n 


dAHs 




(3.18) 


arN'' 




M 


'Mi 


n=0 r{j,n) CTn r=l 


IA(t,r,(Tn) 


i=l 


Fj, 

•j 


(3.19) 

Remark 2 The semigroup property has been used as essential ingredient to simplify the 
above proof. Indeed, even without arguing on the time derivative, we can still restrict to 
such a small interval of time, that only a single collision takes place, in both the A^—particle 
dynamics and the dynamics of proper subsets of particles. This allows in turn to compute 
just the value of 0 —collision and 1 —collision trees, showing compensations among terms 
of type To and terms of type 7j~. In contrast, a control of the full expansion in the 
generic time interval would require to detect more complicated cancellations. Some of 
them have been classihed in [16]. Our proofs, however, do not characterize the complete 
set of compensations leading from the virtual trajectories of the expansion to the unique, 
real motion of particles appearing on the left hand side. 

We provide one example of such compensations, in the case of iV = 3 hard spheres 


already considered in Section 2.3.1 The collisions therein called ci, 02 , 03,04 are all the 


collisions exhibited by the dynamics in [0,t]; see Figure (2.3.11). Let Ai(s) the empirical 
distribution supported on this trajectory. Then the Dirac mass of particle 1 at time t is 
produced in zi{t) by the term r(l, 2) = (1, 2) with (J 2 = (+, +) in the BBGKY expansion 


(i.e. the tree pictured in (2.3.12)). But the same term produces also a virtual Dirac mass 


in a conhguration zi{t) 7 ^ according to (2.3.13). This Dirac mass is compensated 


by a different term of the expansion. That is, the negative contribution coming from 
r(l,2) = (1,2) with (T 2 = (—,+). To hnd this contribution, just let particle 1 ignore the 
last collision with particle 2 in Figure ( [2.3.11 ). 

We deduce now from Theorem [T] the validity of the BBGKY hierarchy for measures: 


Corollary 1 Given an initial measure on with density IFo G L^(AlAr) and invari¬ 
ant for permutations of particles, let {/j(t)}^i be the family of time-evolved marginals 


as defined in (2.1.4)-(2.1.5). Then, the expansion (2.1.12) holds for any t > 0, almost 


everywhere in A4 n- 
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To prove the corollary, let z^it) be distributed according to W(t) and consider Aj(t) 
as a random probability measure over . By Theorem the family of point processes 


satishes (2.1.12): 


^ 3 nt pti 

^(*) = E a{N — j,n) / dti / dt 2 




dtr, 


n=0 


'0 


- t2) ■ ■ ■ 


(3.20) 


for any z^r G Ai%. Thus, we just need to compute the expectation E of the above relation 
with respect to Wq. Observe that this amounts to apply dz^Wo^z^) in both sides. 

, the symmetry in the exchange of particles and ( 3.9[ ), it easily follows 
= fj (t) (p ■ On the other hand, in the assumptions of the corollary. 


Using (|2.1.4|)-(|2.1.5 
that E 




M 


' A^ivxA(t,r,cr„) 


dz]^d(^j^fi 


(Ci+n) 0(Cf(t))lUo(ziv) 


Z=1 


,E 

^1 T” i'^j + n 


f _ ni=i ^is) 

^ iV(iV - 1)... (iV - j - n + 


M 


TT E 0(Cf (f))w^o(zv) 

' i=l 


j+n 


M 


' MNXA{t,r,Crn) 


dzj\fd(^j^n n S{C - Zs) ^(Cf’\t))Wo{zN) 


S=1 


2 = 1 


M 


j+n 


= £ 


^ ^Cj+n dZjj_n n ^(C - Zs) /o, 

A(t,r,(T7j,) 2 = 1 ^ 3^j-\-n S = 1 

« M 

I d(^j^n fo,j+n{Cj+n) 

A{t,r,<Tn) 

^ dQ dA 


,j+n{^j+n) 


where in the last step we applied the change of variables (2.3.2). Therefore, the expectation 


of the right hand side of (3.20) represented in the weak formulation (2.3.3) is equal to the 


(integrated) expansion (2.1.12), i.e. (2.3.4) holds for any bounded continuous 0. This 


concludes the proof of the corollary. 


3.1 Microscopic series solution to the Enskog hierarchy 

In this Section we consider, at a formal level, the Enskog hierarchy and its series solution, 
denoted by for initial microscopic states /iAr(O) = with Ne^ = supported in 
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ztv € which reads (see (2.1.17), (2.2.16), (2.3.6)): 


/i 


ftl rtn-l 

dt2 • • • / dtr, 
'0 Jo 


■(() = ^ A- [ dh [ 

n >0 “^0 “^0 

■Si{t — ti)C2S2{tl — ^ 2 ) ■ ■ ■ Ci+n<Sl+n{tn) 

00 „ 

= E E E^'" / . (3.1.1) 

-n/i ^ ^ 


77=0 r(l,n) CTn 


or also: 


00 


T VTTa.A 


z^ii' 

n=0 r(l,n) (Tn r=l 


—n ^—2n 


lAE(t,r,CTn) 


dCtnf^N 


®(l+r 


’(Cf+J.#'«f(i)). (3.1.2) 


The aim is to understand whether one may conclude that = /iAr(t). Since the 


expansion for satishes the semigroup property (see (2.1.25)), we are allowed to 


proceed exactly as done in the previous section for the H-S hierarchy and focus on the 


single time interval in the partition (3.10). 


In doing so, we shall compare the result (3.19) with the right hand side of (3.1.2). Since 


we fix = A“^, the only difference between the two is the use of the IFF instead of the 
EBF and, correspondingly, the different domains of integration (remind also the discussion 
at the end of Section 2.2). In particular, while the domains in (3.1.2) cover all 


there holds in contrast A{t, F, cr„) C which takes into account the non factorization 


of the initial state in (3.19). Of course, outside Aij, the tensor product 7 ^ 0 because 
of the contractions of type 

5{Ci-z,)5{C2-zJ) . (3.1.3) 


In other words, (3.5) holds only over the phase space “with holes” 


Observe now that the expansions in (3.19) for j = 1 and in (3.1.2) appear to be identical 


when t is so small that Properties 1 and 2 are satished in the single time interval (0, f)[^ 
Indeed in this case there are no recollisions, thus the interacting flow and the Enskog flow 
are the same. However, strictly speaking, the identity fj,^(t) = UNit) remains doubtful 


because the contractions in (3.1.2) give contributions which are absent in (3.19). When 


considering such contributions, we face an ambiguity, as explained by the example that 
follows. 


^The main statement of [1] is in fact that the first equation of the BBGKY hierarchy for empirical 
marginals reduces to the Boltzmann-Enskog equation. But this relies on the interpretation of the second 


marginal at the boundary via the product formula (3.5). 
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Let us consider the simplest case N = 2, Ne^ = A ^ and assume that the two particles 
with initial conhguration Z 2 collide in the time interval (0,t). 


The term n = 0 in (3.1.2) is the free flow As before, this term cancels 

with the negative part of the term n = 1, i.e. which is the contribution due 

to the trees 




Therefore we can write 


— Tq^{C l, + Ti^iCl,■ 


But there are other terms in (3.1.2) which must be evaluated. 
For instance, consider the tree 



(3.1.4) 


The associated contribution, integrated against 0, gives 

/ dC^ dC^ dCi fx%%C^,C2Xi)(l){Ci{'t)) (3.1.5) 

Ja^ 

= Ja^ dCi dCi 0(Cf (t)) 

■ 'd(Cf - Zl) d(C2^ - Z2) diCi^ - Z2) + HCf - Z2) HC2^ “ ^l) ^(Cs^ " ;^l)] , 

where is the image of the EBF-map Cf, t 2 , n; 2 , vi 2 —C3'(0), namely the set defined 
by the condition: “the free flow leads first particle 3 and then particle 2 at e—distance 
from particle 1”. Note that all the other contributions arising from the dehnition of /i®^ 
do vanish identically. 
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The dynamical content of (3.1.5), in terms of the EBF, is the following. Particle 1 
creates particles 2 and particle 3 simultaneously, at times ti = ^ 2 - Hence the deltas are 
evalnated on the border of A^. One can be tempted to interpret the nnmerical valne of 
the above term as just —/iAr(t). This would cancel the contribution (3.1.4). On the other 
hand, by a similar computation, the value of fiNit) would be restored e.g. by 



Furthermore, observe that the term 



corresponds to a triple collision in the Fnskog flow for which Ci{t) is not clearly dehned. 

The ambiguity we are dealing with cannot be easily solved by regularizing the initial 
^-functions. The obtained value of trees would depend on the regularization we choose. 
For instance, a symmetric regularization yields —from (3.1.5), as a consequence 
of the time ordering present in the dehnition of A^. Moreover, families of trees yielding 
the same contribution with alternate signs can be constructed, as in the example above, 
so that the resulting series is divergent (or at most converging in the Cesaro sense). In 
conclusion, it is difficult to give a dehnite meaning to (3.1.2) without further prescriptions. 

It is not possible, however, to solve the problem by giving a prescription that simply 
eliminates contraction^ These are indeed essential to reconstruct the H-S dynamics. In 
fact remind that, in the EBF, the only interactions are creations of particles (see Sec. 


2.2). Now consider for instance a case N = 3 with three or more interactions in the H-S 


dynamics, as in ( |2.3.11 ). In order to generate such a dynamics in (3.1.2), we strictly need 
a term r{j,n) with n > 2, which will have total number of particles j -|- n > 3. Since 


^Or by considering the Enskog hierarchy (2.3.61 with correlated initial data goj = 6j<N2^j- 
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N = 3, this means that at least 2 particles occupy the same configuration. In other words, 
recollisions in the H-S dynamics are replaced by creations and contractions in the EBF. 

One way to eliminate the pathologies we have discussed for the simple example of a 
single collision, is to consider a regularization based on the separation of times in the 
Duhamel series: 


n>0 



■Slit — fi)C2i52(tl — 




(3.1.6) 


With this definition, all the contraction terms discussed above are avoided. In particular, 
the same computation of the proof of Theorem for the case j = 1 shows that, in the 
single time interval of the partition (3.10), = fiN(t) holds. However, the regularized 

series does not converge to /iAr(t) for arbitrary times (notice that the semigroup property 
fails). 


We conclude with the remark that many singular solutions to the Boltzmann-Enskog 
equation which are not corresponding to the particle dynamics may be easily constructed. 
For instance, consider N = 2 and the initial condition 


2 L 


6(5 - xi)S(ti - vi) + 6(5 - x2)S(7]) 


(3.1.7) 


with support outside AI 2 , namely |xi — ;r 2 | < £ ■ Then the unique solution to the equation 

is 


2 L 


6(5 -xi- vit)6(7] - ui) + 6(5 - 12 ) 6 ( 1 ;) 


(3.1.8) 


because only the term n = 0 is non-vanishing. Clearly, many other examples may be 
provided, but it seems hard to construct a solution reproducing the physical dynamics, at 
least in the sense of (2.3.6). 


Acknowledgments. We are grateful to Maxime Hauray for interesting discussions on the 
subject. 
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